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We explain how fluctuations of ratios can constrain and falsify the statistical model of particle 
production in heavy ion collisions, using K/n fluctuations as an example. We define an observ- 
able capable of determining which statistical model, if any, governs freeze-out in ultrarelativistic 
heavy ion collisions. We calculate this observable for K/ % fluctuations, and show that it should 
be the same for RHIC and LHC energies, as well as independent of centrality, if the Grand- 
Canonical statistical model is an appropriate description and chemical equilibrium applies. We 
describe variations of this scaling for deviations from this scenario, such as light quark chemi- 
cal non-equilibrium, strange quark over-saturation and local conservation (canonical ensemble) 
for strange quarks. We also introduce a similar observable capable, together with the published 
K* /K measurement, of ascertaining if an interacting hadron gas phase governs the system be- 
tween thermal and chemical freeze-out, and of ascertaining its duration and impact on hadronic 
chemistry 
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1. Introduction 

One of the main objectives of iieavy ion piiysics is to study the collective properties of strongly 
interacting matter. It's equation of state, transport coefficients, degree of equilibration and phase 
structure, and the dependence of these on energy and system size. 

Thus, the natural approach to study soft particle production in heavy ion collisions is through 
statistical mechanics techniques. Such an approach has a long and illustrious history [|T], ^ ^, ^. 
However, the systematic quantitative comparison of data to the statistical model is a comparatively 
recent field ; A consensus has developed that the statistical hadronization model can indeed fit most 
or all particles for AGS,SPS and RHIC energies ||, g ^ |||[ [T|, [HI]. 



This consensus, is, however, somewhat superficial. While it is true that one can get a reason- 
ably nice-looking fit with a statistical model ansatz, it does not follow that the system is actually 
thermally and chemically equilibrated at freeze-out: Considering the paucity of data points when 
particle abundances are modeled, such a fit is by itself not a guarantee of the physical significance 
of parameters such as temperature and chemical potential. When statistical significance of these 
fits is calculated, it is apparent that the statistical model is nowhere near "proven" according to the 
standards generally accepted in particle physics p2| ] To compound this point, it should be remem- 
bered that a roughly "thermal" distribution is also observed in systems that appear too small to be 



equilibrated, such as p — p and even e^e^ [13]. 

Furthermore, an unambiguous link between heavy ion phenomenology and statistical model 
parameters is still missing. It is unclear weather strangeness enhancement, or features such as the 
horn, are due to canonical effects of incomplete chemical equilibration. It is unclear if |^, and 
at what energy, does light and/or strange quark chemical equilibration occur. It is unclear what 
bias, if any, does in-medium mass modification of short-lived resonances introduce into statistical 
model parameters. A more stringent test of the statistical model could serve both as a strong 
confirmation that statistical physics is an appropriate description for heavy ion collisions, and as a 
tool for deciding which statistical model is more appropriate. 

Particle yield fluctuations are a promising observable to falsify the statistical model and to 
constrain its parameters (choice of ensemble, strangeness/light quark chemical equilibrium) [pl|]. 
One can immediately see that fluctuations are a stringent statistical model test by considering the 
fluctuation of a ratio between two random variables. 

_ {{m?) , ((AA^2)^) , (AyViAyV2) 



Since, for an equilibrated system, ((AA/^i)^) {N) ~ {V), where {V) is the system volume [15], it 
is clear that a^^ depends on the hadronization volume in a manner opposite to that of particle 
yields, inversely rather than directly linearly proportionally. Volume fluctuations (which make 
a comparison between of statistical model calculations to experimental data problematic), both 
resulting from dynamics and from fluctuations in collision geometry, should not alter this very 
basic result since volume cancels out event by event [|l^], provided hadronization volume is the 
same for all particles (a basic statistical model requirement). 

Thus, observables such as (A'^i) o^^^^^, provided (A'^1,2) and o'^^jj^^ are measured using the same 
kinematic cuts, should be strictly independent of multiplicity and centrality, as long as the statistical 
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model holds and the physically appropriate ensemble is Grand Canonical. If the temperature and 
chemical potentials between two energy regimes are approximately the same at freezeout (this 
should be the case for RHIC top energies and LHC, provided chemical equilibrium holds), this 
observable should also be identical across energy regimes. This could be used as a stringent test of 
the statistical model. 

Fluctuations are more sensitive to acceptance cuts than yields. A partial "fix" for acceptance 
cuts that does not require detector-specific analysis is mixed event subtraction, based on the idea 
that fluctuation effects resulting from acceptance cuts are present both in real and mixed events 
(this is the case for fluctuations, but not for correlations). Thus, an appropriate observable to model 



would be [17] 



<^dyn = - (1-2) 

where a,^,-^ is the mixed event width. In the absence of any correlations, it reduces itself to the 
Poisson expectation, 

.2 ^ 11 



i'^mix)Ni/N2 — TTTT + 



putting everything together, the right observable to model becomes equivalent to vf^" 



{Ni) {N2) 

N1/N2 



^dyn _ , dyn \2 _ 2 _ ( PoissonVl _ 
^Ni/N2 ~ ^"Ni/N2> ~^Ni/N2 \^Ni/N2 ' ~ 

(M(M-l)) (jV2(A^2-l)) , (A^lA^2) 



(Vdyn is theoretically equivalent to G^y^, but experimentally it is measured by histogramming []17[]). 



This observable is currently subject of intense experimental investigation [ ]18| , |19| ]. 
We therefore propose to use the scaling of 

= (A^Ovjr^ (1.4) 

to test the statistical model validity among different energy, system size and centrality regimes. 
"Primordial" fluctuations of each observable, ((AA^i 2)^), are calculable from Textbook meth- 



ods [ |15| ] from Fermi-Dirac or Bose Einstein statistics. Observed fluctuations, however, must also 



include corrections from resonance decays. This correction, ((AA/^,-^,)^), is given by 



{{ANj^i)^) = Bj^til -Bj^t) {Nj) +Bj^,{{ANjf) . (1.5) 

where the first term represents the "anti-correlation" due to exclusive decay channels (if the p 
decays into n^n^, it will not decay to 7r"7r*') while the second term is the fluctuation in the number 
of resonances itself. Similarly, the correlation term in Eq. |T]T]is given by 

{ANim)=Y.Bj^n{Nj) (1.6) 

For the analysis described in this section, it is safer to ignore this quantity by choosing particles 
(such as K and 7i~) least correlated by resonances. In the next section we will detail with the 
correlation term as an observable, and show that it is also very useful for distinguishing between 
different freeze-out scenarios. 



3 



Fluctuations 



Giorgio Torrieri 



SHAREv2.X ^ provides the possibility of calculating all ingredients of for any 

hadrons, incorporating the effect of all resonance decays, as well as chemical (non)equilibrium. 

It is important to underline that the value of (calculated from statistical model param- 

eters and, as we will see, sensitive to the degree of chemical equilibration of the system) should 
be constant across any system where the intensive parameters are the same, for example different 
centrality regimes or system sizes at the same energy. For instance, since the chemical potential of 
Cu-Cu 200 GeV collisions should be comparable the chemical potential at Au-Au, 



^K-/7l- 



^ „ — ^K- 71- 



Au—Au 



or in other words 

,dyn 



Cu-Cu '\^k)\au-Au ^ ^'^ 



Au-Au ^ 1^ 



All— All 



(1.7) 



(1.8) 



A large deviation from this value, or a systematic variation of *Fj^J with centrality, should be 
taken as indication of an admixture of non-statistical behavior (for example, a significant effect of 



a non-thermalized "corona" [^, ^). In the presence of a description for the corona, the formulae 
for *F can be easily extended with 

k^)^k^)core^k^) corona (1-9) 
((AyV)2) ^ ((A \}^core+^corona\Y) ^ {{t^coref) + ((AA^™)'> (1-10) 

A break of the flat scaling of would strongly motivate such a description, and modeling it 

could help in constraining the dynamics of the corona. 

The calculation for j.^- , as well as is shown in Fig. |l[ These species were chosen 

because their correlations (from resonance decays, N* NiNj), which would need corrections for 
limited experimental acceptance, are small. 

The methods described in this section can be used to ascertain the origin of the "horn", either 
at NA49 or in the future low energy RHIC runs. If the "horn" is due to a smooth transition between 
baryon-dominated and meson-dominated freeze-out, than should not significantly change 

between energies to the left and the right of the horn's tip. (Hb will change, smoothly, but K and n 
are only relatively sensitive to this). 

On the other hand, if, as hypothesized in |^, the "horn" is due to a sudden jump in the light 
quark phase space occupancy, than n fluctuations should be correspondingly enhanced, resulting 
in a jump of ^^-/j^-- When performing this scan, care should be taken that the yield and the 
fluctuation are measured within the same kinematic cuts (this is the reason why such an analysis is 
not possible with currently available data, but could become possible in future SPS and low energy 
RHIC runs). 

A more quantitative, and striking, signature for deviation from equilibrium can be made be- 
tween RHIC and LHC energies. Equilibrium thermal and chemical parameters are very similar at 
RHIC and the LHC(the baryo-chemical potential will be lower at the LHC, but it is low enough at 
RHIC that the difference should not significantly affect n and K abundance). Thus, ^'^-^j,-- should 
be identical, to within experimental error, for both the LHC and RHIC, over all multiplicities were 
the statistical model is thought to apply. 
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According to |24], chemical conditions at freeze-out (at SPS, RHIC and LHC) deviate from 
equilibrium, and reflect the higher entropy content and strangeness per entropy of the early de- 
confined phase through an over-saturated phase space occupancy for the light and strange quarks 
(Ys > Yij> 1)- If this is true, than ^I'j^J^yy^ should still be independent of centrality for a given energy 
range, but should go markedly up for the LHC from RHIC, because of the increase in 7^ and 7. 

We have calculated ^^'^^j^^ for RHIC and LHC energies, for the sets of parameters used in 
[24|. The left and right panel in Fig. [l] shows what effect three different sets of 7^ ^ inferred in 



[24] would have on ^^-^j^- and In the left panel we have also included the value of 



^K-/7i- ^'^^ '■°P ^i^^i'gy RHIC. As shown in [^], this value for top centrality matches expectations 
for non-equilibrium freeze-out (and is significantly above equilibrium freeze-out). A centrality 
scan of necessary to confirm the consistency of this result has not, however, as yet been 

performed. 

If non-statistical processes (mini-jets, string breaking etc.) dominate event-by-event physics. 



the flat ,^ scaling on centrality /multiplicity should be broken, and ,^ would exhibit a 
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non-trivial dependence on Npan or dN/dy. 
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Figure 1: (color online) ^J-^^- (Left panel) and 1'^-/^+ (right panel) calculated for vaiious statistical 
hadronization parameters [M] at the LHC. The left panel also shows the RHIC calculation | 



This is also true if global correlations persist, such as is the case if the Canonical and micro- 



canonical ensembles [gg, |27|] are physically more appropriate to describe the system than the 
Grand-Canonical ensemble. Whether this is in fact the case is not immediately apparent: Higher 
energy RHIC and LHC experiments are only capable of observing the mid-rapidity region, a small 
fraction of the system where the energy density and chemical potential appear to be constant. The 4 
n acceptance limit (explored at the SPS), where conservation laws have to be obeyed exactly, could 
also correspond to a highly inhomogeneous system (E.g. Baryo-chemical potential is thought to 
vary markedly with rapidity), where no ensemble is appropriate. 

Thus, it would appear that the Grand Canonical ensemble is more appropriate. On the other 



hand, micro-canonical calculations of some SPS observables are remarkably successful [28], sug 
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gesting the possibility that correlations from conservation laws should be strictly taken into accoun 
even within systems where their applicability is not intuitively clear. Furthermore, models have 
appeared in the literature were strangeness is produced and strictly conserved locally, necessitating 
a Canonical approach [|9|]. This theoretical ambiguity makes further experimental study desirable, 
and, due to their ensemble-specificity even in the thermodynamic limit [ p6| , p7| ], fluctuations are a 
very convenient probe. 

If global correlations persist for particle N2 and/or A'^i, than ^'^'^^j^^ becomes reduced, and 
starts strongly varying with centrality in lower multiplicity events. Thus, if strangeness at RHIC/the 
LHC is created and maintained locally, ^'^^ji^^ should develop a "wiggle" at low centrality, and be 
considerably lower than Grand Canonical expectation. For it should be lower by a factor 

of two. 

In conclusion, measuring ^^-^j^- and ^^+^](- , at comparing the results between the LHC and 
RHIC can provide an invaluable falsification of the statistical model, as well as constraints as to 
which statistical model applies in these regimes. 



2. Fluctuations and resonances 



A still unresolved ambiguity of statistical models, with profound repercussions within other 
branches of heavy ion physics is the duration, and impact on hadronic observables, of the phase 
between hadronization (the moment at which particles become the effective degrees of freedom) 
and freeze-out (the moment at which particles stop interacting). 

If chemical freezeout temperature is Tchem ~ 170 MeV, as deduced from equilibrium statistical 
model fits, and thermal freezeout temperature is Ttherm ~ 100 MeV, as deduced from fits to particle 
spectra, it follows that there is a significant interacting hadron gas phase that has the potential of 
altering all soft hadronic signatures. 



The failure to solve the HBT problem [30], combined with acceptable fits obtained by simul- 



taneous freeze-out models [ |31[ , |32| ], suggests however that we are missing something fundamental, 
and more direct probes of freeze-out dynamics should be needed. 



The measurement of Resonance yields offers such a probe |B3, 34], since short-Uved hadronic 



resonances decay before the interacting hadron gas phase (if it exists) is over Thus, rescattering of 
decay products can deplete the amount of observable resonances, while regeneration could create 
additional resonances not present at hadronization. 

], at abundances below equilibrium statisti- 



The observation A(1520) and K*{S92) [||, | 
cal model expectations, could be interpreted as an indication of such reinteraction, with rescattering 
predictably dominating over regeneration. This interpretation, however, is not unique: Chemi- 
cal non-equilibrium fits recover the resonance abundance exactly, with no need for an interacting 
hadron gas phase 



3, 34] 



Thus detection of short lived resonances can not tell us freezeout dynamics unless a different 
signature, more sensitive to chemical freeze-out, is obtained. As is apparent [E 



16] from Eqs 



1.1| and [L6| the correlation term is precisely the required observable, since correlations between 
multiplicities are fixed at chemical freeze-out. This correlation term can be measured by comparing 
observables such as jj^- (correlated by K*{^92)) with jj^- (not correlated by resonances). 
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Figure 2: (color online) ^^j^ and K* /K calculated within the equilibrium and non-equilibrium statisti- 
cal models 
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this result is somewhat spoiled by finite baryochemical potential, as well as higher lying resonances 
(including the anti-correlation term in Eq. [T!^ ). To ascertain the size of these corrections, we have 
used SHARE to calculate both the K* /K^ and A^Fj^^^^. As shown in Fig. | (left panel), these 
corrections make up a 10 % effect, less than the expected experimental error and not enough to 
alter the difference between a single freeze-out and two simultaneous ones. 

A long reinteracting hadron gas phase would in general bring the observed (final) abun- 
dance of K* /K away from the chemical freezeout value (either up, by regeneration, or down, 
by rescattering). Thus, /^^^j^ would become different from K*/K^ (Fig. fright panel). In 
the weak interaction limit, regeneration would presumably be rarer than rescattering so the ob- 
servable K*/K^ abundance would be suppressed by a factor that combines the interaction width 
^regeneration with the lifetime of the interacting phase z 



Ar — ^rescattering 



K* 



observed 



(2.3) 



In the strong reinteraction limit, rescattering and regeneration would reach detailed balance until a 
lower freeze-out Ttherm, so the observed K* /K would be sensitive to the difference between the 
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two temperatures as well as the mass difference (Am) between K* and K 

Am Am 

One untested effect that could spoil this result is strong rescattering capable of bringing the K and 



K* 



observed 



(2.4) 



TT out of the detector's acceptance region in phase space [14, 16]. This effect can not be taken into 
account by mixed event techniques described in the previous section (since mixed events retain no 
two-particle correlations), and calculating it in a model-independent way is problematic. 

Inferring the presence of such a correction is however relatively straight-forward: The proba- 
bility of such rescattering strongly depends on the width of the acceptance region. Thus, if AW^^^ 
as a function of the rapidity window should go from zero (at small rapidity no multiplicity cor- 
relation survives) and saturate at a constant value (where the full resonance derived correlation is 
recovered). This constant value, as long the rapidity window is much smaller than the total extent 
of rapidity of the system, is the quantity that can be related to §^\^fi^^- 

The dependence on centrality of A^'^^^, on the other hand, has to remain flat, since in the 
Grand Canonical limit the ratio of two particles should be independent of centrality, and the total 
system size should not alter the probability of a local process (such as scattering in/out of the 
acceptance region) to occur. 

If A*F^^^ obeys the consistency checks elucidated here (saturating value when rapidity window 
is varied, flat when centrality is varied), it should be taken as a reliable indication that the A^'^^^ 
measurement in fact reflects the value of ^ at chemical freezeout. 

If A*F^y^ depends on rapidity up to the acceptance region of the detector, a more thorough 
effort to account for acceptance corrections to the correlation is needed. This can be done by using 



the same techniques utilised for direct resonance reconstruction [|35|, However, such an 

endeavour is detector specific, and goes well beyond this write-up. 

In conclusion, we have shown that observables incorporating both yields and fluctuations give 
a stringent test of statistical models. We have also argued that such observables can be used to 
infer the duration of the interacting hadron gas phase, and its effect on hadronic observables. We 
expect that forthcoming experimental data, together with the methods elucidated here, will allow 
us to clarify some of the outstanding puzzles apparent in the study of heavy ion collisions. 

GT thanks the Alexander Von Humboldt foundation, the Frankfurt Institute for Theoretical 
Physics and FIAS for continued support We would also like to thank Gary Westfall, Sangyong Jeon, 
Marek Gazdzicki, Mike Hauer, Johann Rafelski and Mark Gorenstein for useful and productive 
discussions. 
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